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ABSTRACT 


Static pressure distributions, mean velocity profiles, 
distributions of turbulence intensities, and Reynolds 
stresses, spatial correlation coefficients, and the Taylor 
microscale of turbulence were measured across the stern 
boundary layer of an axisymmetric body with an inflected 
stern. A theoretical and numerical iteration scheme, 
which uses the boundary layer and open wake displacement 
body, is found to model satisfactorily the viscid-inviscid 
interaction between the thick stern boundary layer and the 
external potential flow. The measured values of turbulence 
intensity, eddy viscosity, and mixing-length parameters in 
the stern region are found to be much smaller than those of 
a thin boundary layer. An approximate similarity character- 
istic for the thick axisymmetric stern boundary layer is 
obtained when the measured mixing-length parameters, the 
measured correlation length scales, and the measured Taylor 
microscales are normalized by the square root of the 
boundary-layer cross-sectional area instead of the local 
boundary layer thickness. When this simple similarity 
hypothesis for the mixing length and the displacement body 
is incorporated into the McDonnell Douglas Corp., Cebeci- 
Smith differential boundary-layer method, modified to 
consider the displacement body and wake, the theory pre- 
dicts very well the measured distributions of the mean 
velocity throughout the entire stern boundary layer. The 
computation method developed earlier was found to predict 
well the boundary layer on convex sterns and is now found 
to apply equally well to the flow on inflected sterns. 


ADMINISTRATIVE INFORMATION 
The work described in this report was funded under the David W. Taylor 
Naval Ship Research and Development Center's Independent Research Program, 


Program Element 61152N, Project Number ZR 000 01, and Work Unit 1552-103. 


INTRODUCTION 
This report presents measurements of the turbulent boundary layer 
characteristics in the thick stern boundary layer of an axisymmetric body 
with an extremely full inflected stern. These measurements include the 
distributions of pressure coefficient and shear stress on the body, static 
pressure distribution across the stern boundary layer, mean velocity pro- 


files, turbulence intensities, Reynolds stresses, eddy viscosity, mixing 


length, and spatial correction function. This investigation is an ex- 
tension of the earlier research on convex axisymmetric thick stern boundary 
layers reported by Huang, et aires 

The experimental measurements are compared with the results computed 
by the method of Wang and Boner for computing the potential flow- 
boundary layer interaction on axisymmetric bodies. The initial step in 
this procedure begins with the computation of the potential-flow pressure 
distribution on the body using the computer code of Hess and Smet.” The 
McDonnell Douglas Corp., Cebeci-Smith** methods is then used to calculate 
the flow over the body and the integral relations of @ranvaliies are used to 
calculate the flow in the wake. The local displacement thickness computed 
in these boundary layer calculations is added to the original body to 
obtain a modified body and wake geometry as suggested by Peeecend and 
De neneells The procedure is repeated until the pressure distributions on 
the body from two successive approximations agree to within a given error 
criterion. This iterative scheme is referred to as the displacement body 
concept. 

Several modifications have been made to the Douglas C-S differential 
boundary layer meenode by Wang and Huenenes These changes were incorporated 
to improve the comparison of the thick boundary layer parameters with the 
experimental data of Huang et mi and Patel and teas” Transverse 
curvature effects are included in the solution of the standard thin 
boundary layer equations in the Wang and Huang method. Effects due to 
longitudinal curvature are neglected. The pressure variation across the 
boundary layer is obtained from potential flow calculations for the final 
displacement body. The mixing length and eddy viscosity in the thick 
boundary layer region are obtained by modifying the thin boundary layer 
values by factors which are proportional to the annular area between the 
body surface and the edge of the boundary layer. 

In the stern and near-wake region (0.95 < x/L eO5))5 where x is the 


axial distance from the nose and L is the total body length, a fifth degree 


*A complete listing of references is given on page 51. 


**The McDonnell Douglas Corp., Cebeci-Smith method will hereafter be 
designated as Douglas C-S. 


polynomial was used to determine the boundary layer flow. The coefficients 
for this polynomial were determined by the conditions that the thickness, 
slope, and curvature be equal to those calculated by the modified boundary- 
layer method at x/L = 0.95 and by the Granville integral wake wallatsones® at 
x/L = 1.05. This method was found by Huang et ail,” to give excellent 
agreement with experimental values of pressure, shear stress, and velocity 
profiles over the forward 90 percent of the bodies investigated. As the 
boundary-layer thicknesses became larger, particularly when these thick- 
nesses became greater than the radii of the bodies, the measured values of 
shear stress and velocity became smaller than those predicted by the 
theory. 

In the following discussion, the experimental techniques and 
geometries of the model are given in detail. The experimental and theo- 
retical results are compared and the measured turbulence characteristics 
are presented. The method of obtaining the eddy viscosity and mixing 
length is discussed and the experimentally—determined distributions are 
given. The concept of similarity length-scale in the thick stern boundary 
layer is examined experimentally. The square root of the annular area be- 
tween the body surface and the edge of the boundary layer is found to be 
the appropriate length scale for the axisymmetric thick stern boundary 


layer. 


WIND TUNNEL AND MODEL 
The experimental investigation was conducted in the DINSRDC anechoic 
wind tunnel facility. The wind tunnel has a closed jet test section that 
is 8 ft (2.4 m) square and 13.75 ft (4.19 m) long. The corners have 
fillets which are carried through the contraction. The test section is 
followed by an acoustically-lined large chamber 23.5 ft (7.16 m) long. It 


was found previously by Huang et Ate that the ambient free-stream 


turbulence levels, (/ u'~/U,) x 100, are 0.075, 0.090, 0.100, and from 
0.12 to 0.15 for free-stream velocities U5 Oe ZHGAL S0555 Sk3odl, Eioxsl “AS 7/ 
m/s, respectively. Integration of the measured noise spectrum levels in 


the test section from 10 to 10,000 Hz indicated that the typical background 


acoustic noise levels at 30.5 m/s were about 93 dB re 0.0002 dyaviens 
(0.0002 Pa). These levels of ambient turbulence and acoustic noise were 
considered low enough so as not to unfavorably affect the measurements of 
boundary-layer characteristics. The maximum air speed which can be 
achieved is 200 ft/sec (61 m/s); in the present experiments the wind tunnel 
velocity was held constant at 150 ft/sec (45.72 m/s). 

An axisymmetric inflected (convex and concave) afterbody with a bow- 
entrance length-to-diameter ratio (L/D) of 2.05 was used for the present 
experimental investigation. The detailed offsets are given in Table 1. 
Figure 1 shows a comparison of the present afterbody (Afterbody 5) and 
Afterbodies 1 and 2 previously investigated by Huang et ailer 2? The contour 
of Afterbody 5 changes from convex to concave at x/L = 0.913, whereas 
Afterbodies 1 and 2 are convex up to x/L = 0.96. Afterbody 5 is connected 
to a parallel middle body which is 6.0 ft (1.83 m) long and an existing 
streamlined forebody with a bow-entrance length-to-diameter ratio (L/D) of 
1.82. The total model length is 9.55 ft (2.91 m) with a maximum forebody 
diameter of 0.917 ft (0.28 m). 

The model was supported by two streamlined struts separated by one- 
third of the model length. The upstream strut had a 15-cm chord and the 
downstream strut a 3-cm chord. The disturbances generated by the support- 
ing struts were within the region below the horizontal centerplane, there- 
fore, all of the experimental data were taken above the model on the 
vertical centerplane along the upper meridian where there was little effect 
from the supporting struts. One-half of the model length protruded beyond 
the closed-jet working section into the open-jet section. The ambient 
static pressure coefficients across and along the entire open-jet chamber 
(7.2 m X 7.2 m X 6.4 m) were found to vary less than 0.3 percent of the 
dynamic pressure. Tunnel blockage and longitudinal pressure gradient 
effects along the tunnel length were almost completely removed by testing 
the afterbody in the open-jet section. 

The location of the boundary-layer transition from laminar to turbu- 
lent flow was artificially induced by a 0.024 in. (0.61 mm) diameter trip 


wire located at x/L = 0.05. Huang et ails found that the trip wire 


TABLE 1 - MODEL OFFSETS 


x/L ro/L O/T max x/L Tg trey ire 
.4450  .0480 1.0000 .7768  .0480 1.0000 
-4659  .0480 1.0000 .7821 .0480 1.0000 
. 4868 -0480 1.0000 - 7863 -0480 1.0000 
"5078  .0480 1.0000 7905 .0480 1.0000 
15256 .0480 1.0000 .7936  .0480 1.0000 
~5402 -0480 1.0000 7967 -0920 1.0000 
15518  .0480 1.0000 8000 .0480 +.0000 
15601 .0480 1.0000 .8030 .0480 1.0000 
(5675  .0480 1.0000 8062 0480 .9998 
15727 .0480 1.0000 .8093. .0479. .99¢1 
15779 .0480 1.0000 8135 .0479 9975 
Poe o2 -0480 1.0000 8187 0477 - 9937 
.5884 -0480 1.0000 8256 .0473 9852 
'5936 .0480 1.0000 8308  .0468 .9756 
15989. .0480 1.0000 .8360 .0462 .9632 
.6041 -0480 1.0000 8413 -0455 -9478 
.6093 -0480 1.0000 -8465 -04496 9296 
.6146 -0480 1.0000 -8517 -0436 9085 
"6198  .0480 1.0900 .8570 .0425 8847 
-6250  .0480 1.0000 .8622 .0412 +.8583 
-6303  .0480 1.0000 .8674 0398 .8294 
(6355 .0480 4.0000 8727 .0383 .7983 
16407 .0480 1.0000 .8779  .0367 7652 
"6460  .0480 1.0000 .8831  .0350 .7303 
-6512  .0480 1.0000 8884 .0333 .6939 
(6564  .0480 1.0000 8936. .¢315 .6561 
16617 .0480 1.0000 .8988 .0296 6174 
"6669 .0480 1.0000 9041 .0277. +.5778 
"6722 .0480 1.0000 9093 .0258 .5378 
16774  .0480 1.0000 .9145  .0239 .4975 
"6826  .0480 1.0000 9198 .0219 .4574 
16879  .0480 1.0000 .9250 .6200 .4176 
+6931 .0480 1.0000 9303 .0182 .3784 
“6983 .0480 1.0000 .9355 .0163 .3402 
. 7036 -0480 1.0000 -~9407 0145 - 3032 
17088  .0480 1.0000 .9460 .0129 2678 
/7140  .0480 1.0000 9512 .0113.—-. 2345 
17193. .0480 1.0000 9564 .0098 .2025 
17245  .0480 1.0000 9617 .0084 1757 
|7297 .0480 1.0000 9659 .0073. +. 1517 
-7350  .0480 1.0000 9721 0064 .1326 
"7402 .0480 1.0000 .9774  .0058 1199 
17454 .0480 1.0000 .9826 .0055 .1149 
17507 .0480 1.0000 9832 .0055 .1148 
"7859  .0480 1.0000 (9874 .0053.—-. 1108 
"7611 +.0480 1.0000 .9916  .0048 .0994 
17664 .0480 1.0000 .9958 .0036 .0744 
|7716 .0480 1.0000 1.0000 0.0000 0.0000 
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Figure 1 - Three Afterbodies 


effectively tripped the flow 0.39 in. (1 cm) downstream from the wire and, 
because of the finite parasitic drag of the wire, the boundary-layer can be 
considered to become turbulent at a virtual origin 5.4 in. (13.7 cm) up- 
stream from the trip wire at a Reynolds number of 9.3 x 10 The virtual 
origin for the turbulent flow is defined such that the sum of the laminar 
frictional drag from the nose to the trip wire, the parasitic drag of the 
trip wire, and the turbulent frictional drag aft of the trip wire is equal 
to the sum of the laminar frictional drag from the nose to the virtual 
origin and the turbulent frictional drag from the virtual origin to the 
after end of the model (similar to Reference 10). The virtual origin 


location was taken as the transition location at which the boundary-layer 


changed from laminar to turbulent in the boundary-layer calculations. 


INSTRUMENTATION 
A Preston tube with a 0.072 in. (1.83 mm) inside diameter was attached 
and alined with the flow on the upper meridian of the stern to measure the 
shear stress distribution. The Preston tube was calibrated in a 1 in. 
(2.54 cm) diameter water-pipe flow facility described by Huang and von 
Keeney A series of 0.031 in. (0.8 mm) diameter pressure taps were 
embedded on the upper meridian of the stern at the Preston tube locations. 


These pressure taps were connected to a multiple pressure scanivalve system 


that takes one integral pressure transducer with its zeroing circuit and 
measures a single pressure in sequence along the stern upper meridian. 

The pressure transducer used was designed for measuring low pressure up to 
O62 per) CLos7o: x 10° Pa). The zero-drift linearity, scanivalve, hystere- 
sis, and pressure transducer with zeroing circuit were carefully checked 
and the overall accuracy was found to be within 0.5 percent of the dynamic 
pressure. 

A Prandtl-type static pressure probe of 0.123 in. (3.125 mm) diameter 
with four equally spaced holes located three probe diameters after the 
probe nose was used to measure the static pressure across the boundary 
layer. By yawing the probe in the free stream, it was found that the 
measured static pressure was insensitive to the probe angle up to a 5- 
degree yaw. In all cases, the static pressure probe remained in alinement 
by less than 5-degree yaw. 

The mean axial and radial velocities and the turbulence intensities 
for the Reynolds stress calculations were measured by a TSI, Inc. Model 
1241-20 "X" type hot-film probe. The probe elements are 0.002 in. (0.05 
mm) in diameter with a sensing length of 0.04 in. (1.0 mm). The spacing 
between the two cross elements is 0.04 in. (1.0 mm). A typical schematic 
of the hot-film probe used is shown in Figure 2. A two-channel hot-wire 
and hot-film anemomcter with linearizers was used to monitor the response 
of the hot-film probe. A temperature compensating sensor (probe) was used 
with each hot-film element to regulate the operating temperature of the 
sensor with the changes in air temperature. The "X" hot film and its 
temperature-compensated sensor were calibrated together through the ex- 
pected air temperature range and supplied with their individual lineariza- 
tion polynomial coefficients at the factory. A single element TSI, Inc. 
Model 1212-TI.5 hot-wire probe also was used to measure the velocity 
distribution across the stern boundary layer. Two single-element TSI 
(Model 1212) hot-wire probes with a 90-deg bend separated by a small radial 
distance also were used to measure spatial correlation functions in the 


thick boundary layer. 
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Figure 2 - Schematic of a Two Element Sensor Alined 90 Degrees 
to Each Other and 45 Degrees to Probe Axis 


The frequency response of the anemometer system, for reliable measure- 


ments claimed by the manufacturer, is 0 to 200 
hot film was made before and after each set of 
that the hot-film anemometer system had a +0.5 
ft/sec (+0.23 m/s) accuracy at the free stream 
(45.72 m/s), during the entire experiment. An 


cross-flow velocity by yawing the "X" hot-film 


kHz. Calibration of the “X" 
measurements. It was found 
percent accuracy, +0.75 
velocity of 150 ft/sec 
estimate was made of the 


probe in the freestream. It 


was found that the crossflow velocities were about one percent of the 


freestream velocity. 


The linearized signals were fed into a Time/Data Model 1923-C real- 


time analyzer. Both channels of the analog signal were digitized at a 


rate of 128 points per second for 8 seconds. 


These data were immediately 


analyzed by a computer to obtain the individual components of mean 


velocity, turbulence fluctuation, and Reynolds 


stress on a real time basis. 


A traversing system with a streamlined strut was mounted on a guide 
plate that permitted the traverse to be locked in various stationary 


positions parallel to the longitudinal model axis. 


COMPARISON OF EXPERIMENTAL AND THEORETICAL RESULTS 
Two coordinate systems for an axisymmetric body are given in Figure 3. 


The coordinate system used to experimentally measure the boundary-layer 


Figure 3 - Axisymmetric Coordinate Systems 


flow and to report the results is the x-r coordinate system. The axial 
coordinate x is measured from the nose of the body parallel to the axis of 
revolution. The radial component r is measured from the axis of revolution 
and is normal to the x-axis. The curvilinear s-n coordinate system is used 
in the Douglas C-S differential boundary layer echoed The arclength 
coordinate s is measured parallel to the body meridian and the tangent 
coordinate n is measured normal to the body meridian. 

The Douglas C-S boundary layer method is used in conjunction with the 
displacement body concept to predict the flow over the body. Wang and 
Beene have made several modifications to the Douglas C-S code in an 
attempt to better model a thick stern boundary layer. A comparison of 
their original theory with experimental results is given by Huang et pil 
That comparison revealed several inadequacies in the theory. Wang and 
Teme” subsequently modified their original method based on these experi- 
mental comparisons. The new modifications are evaluated in the present 
investigation. The mixing length in the thick stern boundary layer region 
is modeled by a simple algebraic formula proposed by Huang et an This 


formulation includes only transverse curvature effects. The eddy viscosity 


in the outer region was also modified by an algebraic formula with the in- 
tention of improving the accuracy of the velocity profiles near the body 

in the thick boundary layer region. The final modification to the original 
method uses the predicted potential flow velocities outside of the dis- 
placement body to modify the velocities in the direction parallel to the 


body meridian tangent. 


MEASURED AND COMPUTED PRESSURE AND SHEAR 
STRESS DISTRIBUTIONS 


Pressure taps were used to measure the steady pressure at 15 selected 
locations along the stern surface of the axisymmetric body. The pressure 


coefficient, on is computed from the measured pressures by the relation- 


ship 
PePy P12) 
2 pep, i - m 
Eras 7 PUG 
where p = measured local static pressure 


Pg = measured ambient pressure 
Dr measured dynamic total pressure 
p_ = measured static pressure 
Oo = mass density of the fluid 
The analytically-predicted pressure coefficient was computed on the dis-— 


placement body using the Douglas-Neumann potential flow method of Hess and 


NZ 
Cpa =) (2) 
P (x 


where Ws is the computed potential flow velocity on the displacement body 


smith’ and is given by 


and Uy is the freestream velocity, 150 ft/sec (45.72 m/s). 


10 


A comparison of the measured and theoretical values of the pressure 
coefficient on the surface of the axisymmetric body under consideration is 
given in Figure 4. The theory underpredicts the magnitude of the minimum 


C value slightly; the theoretical trough is not as full as the measured 
P 
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Figure 4 - Computed and Measured Stern Pressure Distribution on 
Afterbody 5 


trough. The three measured values of Ge for = 0.975 fluctuate somewhat 
about the theoretical curve for Sos These discrepancies are small, and 
overall agreement between theoretical and measured values of the pressure 
coefficient is considered excellent. 

Readings from a Preston tube, which was taped to the stern surface at 
the pressure tap locations, were used in conjunction with the steady 


pressure readings to obtain the shear stress distribution at the body sur- 


face. The calibration curve, presented by Huang and von Kerecelan fora 


IL 


Preston tube in a pipe was used to compute the shear stress at the wall, 


Cee The shear stress coefficient Ge is given by 


2 
qc U 
Gi yee (3) 
ESS es 


for the measured shear stress and by 


ANS (anne (4) 


for the analytical shear stress. 

Figure 5 shows a comparison of the measured and theoretical values of 
the shear stress coefficient. Agreement is good for values of = <a) ole 
The trough, which occurs ioe = > 0.9, is somewhat overpredicted by the 
theory. This discrepancy in the distribution of the shear stress caq- 
efficients indicates that the present analytical model may be inadequate 
for the precise prediction of the shear stress on the last ten percent of 


the body. The measured values of C and Cr are tabulated in Table 2. 


MEASURED AND COMPUTED STATIC PRESSURE 
DISTRIBUTION 


The measured and computed static pressure coefficients for Afterbody 5 
are compared in Figure 6 at various locations across the stern boundary 
layer. The off-body option of the Douglas potential-flow computer code was 
used to compute the static pressure distributions for the displacement body 
(solid lines) and the actual body (broken lines). As can be seen in | 
Figure 6, except at x/L = 0.704, 0.987, and 1.045, the measured static 
pressure distributions agree better with the theoretical pressure distri- 
butions computed from the displacement-body model than those from the 
original body. At x/L = 0.704, 0.987, and 1.045, the two calculation 


methods agree to within one percent. The discrepancy between the measured 


eZ, 
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R, =9.3 X 10 
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Figure 5 - Computed and Measured Shear Stress Distribution 
on Afterbody 5 
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TABLE 2 - MEASURED PRESSURE AND SHEAR STRESS COEFFICIENTS 


So Oo CoC © G2 CoO Lo2] FC LSC CS EC LC LoCo SCS 


is 
Iho 
dhe 
ili 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
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Figure 6 - Computed and Measured Static Pressure Distributions Across Stern 
Boundary Layer of Afterbody 5 
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Figure 6a - Nondimensional Axial Lengths, x/L = 0.704, 0.831, 0.873 
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Figure 6 (Continued) 
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Figure 6b - Nondimensional Axial Lengths, x/L = 0.951, 0.909 
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Figure 6c - Nondimensional Axial Lengths, x/L = 0.987, 1.045 
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and computed (by displacement-body model) values of e is, in general, less 
than 0.01 which is about the accuracy of the measurement. The measured 
values of C_ are tabulated in Table 3. 

As gal be seen later, the displacement body of Afterbody 5 is convex 
from the parallel middle body up to x/L = 0.898 and becomes concave down- 
stream from x/L > 0.90, which is rather close to the inflection point of 
the actual body (x/L = 0.913). The differences in the theoretical values 
of off-body pressure distributions computed from the displacement body and 
the actual body are very noticeable at x/L = 0.951 where the measured 
pressure distributions agree rather well with the distributions computed 
from the displacement-body model. The close agreement between the computed 
displacement-body and measured static pressure distributions for Afterbody 
5, as well as Afterbodies 1 and 2 (Reference 1), support the displacement- 
body concept for computing the potential flow outside of the displacement 
surface. The pressure coefficients inside the displacement body may be 
more accurately calculated by a streamline curvature method proposed by 
Deen However, in the present calculation method the value of c. is 
assumed to be constant inside the displacement body. The maximum error of 
the value of S&S from this approximation is less than 2 percent inside the 


displacement body (Figure 6). 


MEASURED AND COMPUTED MEAN VELOCITY PROFILE 

MEASUREMENT 

The axial and radial mean velocity components were measured by the 
"xX" hot-film sensor. At a given axial location, measurements were taken 
at numerous radial positions to obtain the velocity profile. The hot-film 
probe was mounted with the sensing elements in a vertical position. At 
each radial position, 1024 data values were collected in 8 seconds. The 
linearized data were input to a computer code which calculated and printed 
the values of the mean velocity components uy. in the axial direction and 
Mie in the radial direction (given in Table 4). 

A single element hot-wire sensor also was used to compare with the 


results of the "X'" sensor. The values of the axial mean velocity component 
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TABLE 3 -— MEASURED STATIC PRESSURE COEFFICIENTS ACROSS STERN 
BOUNDARY LAYER AND NEAR WAKE 


x/L = 0.7036 x/L = 0.8308 
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TABLE 3 (Continued) 


oe 
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TABLE 3 (Continued) 
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TABLE 3 (Continued) 


x/L = 1.045 
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TABLE 4 (Continued) 


x/L = 0.951 ve ire = 0.2346 tan a = - 0.2952 
Q° max 
T T 
LTE@ mx Vr Yen |) sg |p Wecane u'v' | -u'v! 2 S@ € L 
v, Ti. U U uy | 2 5 
rmax 0 0 0 0 0 U5 q r U56 V ( +6 22 2 
P Tor 0) 
IL 
0.0236 | 0.251 -0.050 0.036 0.026 0.030 0.0443 0.154 0.0347 0.0004 0.0054 0.0041 
0.0447 0.281 -0.053 0.041 0.028 0.032 0.0584 0.167 0.0657 0.0011 0.0127 0.0098 
0.0588 0.297 -0.055 0.044 0.032 0.034 0.0775 0.191 0.0864 0.0015 0.0178 0.0137 
0.0693 0.322 -0.059 0.045 0.032 0.035 0.0852 0.198 0.1018 0.0024 0.0219 0.0169 
0.0869 0.345 -0.061 0.047 0.033 0.037 0.0830 0.178 0.1277 0.0026 0.0245 0.0188 
0.1009 0.374 -0.064 0.048 0.035 0.038 0.0924 0.189 0.1482 0.0026 0.0226 0.0174 
0.1220 0.413 -0.070 0.051 0.036 0.039 0.1053 0.194 0.1792 0.0027 0.0227 0.0174 
0.1431 0.462 -0.077 0.052 0.036 0.041 0.0996 0.177 0.2102 0.0026 0.0223 0.0172 
0.1712 0.511 -0.079 0.050 0.036 0.041 0.1077 0.190 0.2515 0.0029 0.0235 0.0181 
0.2134 0.583 -0.082 0.049 0.036 0.041 0.0973 0.182 0.3135 0.0027 0.0234 0.0180 
0.2275 0.612 -0.085 0.044 0.032 0.041 0.0744 0.161 0.3342 0.0024 0.0236 0.0182 
0.2521 0.638 -0.087 0.044 0.031 0.040 0.0729 0.161 0.3703 0.0035 0.0349 0.0269 
0.2766 0.667 -0.086 0.043 0.032 0.040 0.0791 0.179 0.4063 0.0038 0.0366 0.0282 
0.2977 0.690 -0.085 0.042 0.030 0.039 0.0664 0.160 0.4373 0.0036 0.0369 0.0284 
0.3329 0.723 -0.086 0.040 0.030 0.038 0.0672 0.170 0.4890 0.0041 0.0423 0.0326 
0.3996 0.777 -0.083 0.038 0.029 0.034 0.0568 0.167 0.5870 0.0040 0.0446 0.0343 
0.4734 0.833 -0.080 0.034 0.024 0.027 0.0436 0.175 0.6954 0.0027 0.0346 0.0268 
0.5507 0.887 -0.077 0.023 0.015 0.016 0.0130 0.130 0.8090 0.0014 0.0331 0.0255 
0.6421 0.913 -0.073 0.009 0.007 0.007 0.0001 0.003 0.9432 0.0000 0.0100 0.0078 
0.7475 0.925 -0.062 0.004 0.003 5.003 
0.9233 0.931 -0.053 0.003 0.002 0.003 
1.0322 0.934 -0.047 0.003 0.002 0.002 
1.4539 0.950 -0.030 0.003 0.003 0.002 
1.8123 0.961 -0.022 0.003 0.001 0.002 
2.2024 0.963 -0.012 0.003 0.003 0.002 
2.6662 0.976 -0.002 0.003 0.002 0.002 
x k 
6 6. 6 
—P-_ = 0.1952, = 0.2322, F_ = 0.654 
max max max 
8 nN Us 
—P_ = 0.1000, — >= 0.0480, — = 0.935 
r 2 U 
max r 0 
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TABLE 4 (Continued) 


x/L = 0.987 TE eee = 0.1109 tan a = - 0.0899 
9) Ux | vy | V2 V2 V2 100 -u'v' -u'v' SO € ge L 
* 
ee Uo Yo Uo Uo Uo ua a o, U8, ce V pf6)2e? 
0.0182 | 0.365 | -0.014 | 0.031 | 0.019 | 0.026 | 0.0132 0.065 | 0.0233 | 0.0001 | 0.0024 0.0021 
0.0252 | 0.371 | -0.016 | 0.031 | 0.022 | 0.026 | 0.0196 0.092 | 0.0324 | 0.0001 | 0.0144 0.0127 
0.0322 | 0.378 | -0.018 | 0.032 | 0.025 | 0.026 | 0.0227 0.098 | 0.0414 | 0.0010 | 0.0161 0.0142 
0.0533 | 0.386 | -0.021 | 0.035 | 0.026 | 0.028 | 0.0372 0.140 | 0.0684 | 0.0017 | 0.0206 0.0182 
0.0604 | 0.395 | -0.021 | 0.037 | 0.027 | 0.028 | 0.0410 0.135 | 0.0775 | 0.0019 | 0.0217 0.0191 
0.0814 | 0.412 | -0.025 | 0.041 | 0.030 | 0.031 | 0.0564 0.160 | 0.1045 | 0.0025 | 0.0254 0.0224 
0.1025 | 0.434 | -0.025 | 0.042 | 0.031 | 0.033 | 0.0633 0.168 | 0.1316 | 0.0029 | 0.0269 0.0238 
0.1306 | 0.465 | -0.026 | 0.049 | 0.032 | 0.035 | 0.0757 0.164 | 0.1677 | 0.0034 | 0.0294 0.0259 
0.1482 | 0.487 | -0.028 | 0.049 | 0.034 | 0.036 | 0.0890 0.185 | 0.1903 | 0.0039 | 0.0313 0.0276 
0.1763 | 0.520 | -0.028 | 0.050 | 0.035 | 0.037 | 0.0871 0.171 | 0.2263 | 0.0038 | 0.0307 0.0271 
0.2290 | 0.593 | -0.030 | 0.052 | 0.033 | 0.036 | 0.0828 0.164 | 0.2940 | 0.0036 | 0.0294 0.0260 
0.2853 | 0.655 | -0.031 | 0.046 | 0.030 | 0.036 | 0.0629 0.147 | 0.3662 | 0.0032 | 0.0305 0.0269 
0.3520 | 0.721 | -0.030 | 0.042 | 0.029 | 0.035 | 0.0541 0.141 | 0.4549 | 0.0032 | 0.0325 0.0287 
0.4610 | 0.809 | -0.031 | 0.038 | 0.027 | 6.033 | 0.0524 0.161 | 0.5918 | 0.0042 | 0.0432 0.0381 
0.5383 | 0.855 | -0.029 | 0.035 | 0.025 | 0.029 | 0.0421 0.158 | 0.6910 | 0.0037 | 0.0431 0.0380 
0.6085 | 0.899 | -0.028 | 0.029 | 0.019 | 0.021 | 0.0224 0.135 | 0.7811 | 0.0023 | 0.0370 0.0326 
0.6823 | 0.929 | -0.025 | 0.015 | 0.011 | 0.013 | 0.0056 0.109 | 0.8759 | 0.0011 | 0.0359 0.0317 
0.7456 | 0.938 | -0.025 | 0.006 | 0.005 | 0.009 
0.8124 | 9.941 | -0.022 | 0.004 | 0.004 | 0.006 
0.8897 | 0.946 | -0.020 | 0.003 | 0.003 | 0.003 
1.0583 | 0.955 | -0.017 | 0.006 | 0.003 | 0.003 
1.2305 | 0.959 | -0.013 | 0.005 | 0.003 | 0.003 
1.5890 | 0.961 | -0.006 | 0.004 | 0.004 | 0.003 
1.9369 | 0.971 | -0.001 | 0.004 | 0.003 | 0.003 
2.3094 | 0.977 | 0.001 | 0.004 | 0.003 | 0.004 
2.6397 | 0.977 | 0.008 | 0.004 | 0.003 | 0.004 
2.8013 | 0.979 | 0.010 | 0.004 | 0.004 | 0.004 
on § 8 
aaa = 0.1960, =305257/5 5 F = 0.745 
max max max 
8 U 
—?_ = 0.1103, As 0.0%, “26 aoa 
vmax r a YO 
max 


Zl 


measured by the single-element probe are presented in the figures along 


with the values measured by the two-component probe. 


THEORY 

The Douglas C-S acinad solves the standard incompressible, steady 
continuity and momentum equations for thin axisymmetric boundary layers. 
This set of partial differential equations, with the addition of the 


transverse curvature effects, is 


d(ru.) d(rv_) 


——————_— + —— = 
ds on : (5) 
and 
du du du 
mss ues) Sel. dpe lpm HESS Secrecy 
“s 0s Aol pion p ds Hy Som |S On “s/n (6) 
where u = mean velocity component in the s direction (parallel 
Ss Sai 
to the body meridian) 
Vin 7 mean velocity component in the n direction (normal to 
the body meridian) 
Pp = pressure on the body 
= rQ(s.n) +n cos a 
ro = body radius 
x = axial distance measured from the nose 
Vv = kinematic viscosity of the fluid 
u'v! = Reynolds stress 
sn 
ul = velocity fluctuation in the s direction 
Ue = velocity fluctuation in the n direction 
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The wall boundary conditions are u., (s,n) = va (s,n) = Oat n=0. Effects 


due to longitudinal curvature are neglected. The Reynolds stress model 


is given in detail in the section presenting the turbulence characteristics. 


The Reynolds stress pee in the modified Douglas C-S computer code 


Spates 


is modeled by 
eC. for O<n<n 
il — ce 
du. 
Seite == 
uve wn (7) 
€ ore iy SG i SO 
0 CH= = 
du 
2 NSE Ss a 5 P ; 
where €, = & i) — = eddy viscosity in the 
i r on ; ; 
0 inner region 
* 
EQ = 0.0168 U, o, ieee = eddy viscosity in the 


outer region 


r 
2 = 0.4 to Qn Leal l-exp - =~ in all = mixing length parame- 
ro 0 ter in the inner region 
pa -1/2 
A= 26v (= = Van Driest's damping 
P factor 
a= = i = 
6 89.995 Beret layer thick 
Soe 
BT 0.995 
e 
* : eS 
6. = (1- 3° Jan = planar displacement 
P 0 e thickness 
U = computed potential flow 
e : 
velocity at the edge 
of the displacement 
body 
e = wall shear stress 
n. = value of n at which 
S, = & 
a 0 
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The Douglas C-S computer code is used to calculate the velocity 
components u_ and We ie aiae to and normal to the body surface, respec- 
tively. Wang and Huang modified these velocity components by using the 
predicted potential flow velocities outside of the displacement body and 
converted the normal and tangential velocity components to axial and radial 


components for comparison with the experimental values. 


COMPARISON OF MEAN PROFILES 

Figure 7 shows a comparison of the theoretically predicted and 
experimentally measured mean axial uy and radial We velocity profiles at 
various axial locations on the body. The velocity components are non- 
dimensionalized by the free stream velocity Up: The solid curves represent 
the theoretical results, the open symbols are the ''X" hot-film values, and 
the solid symbols are the single element hot-wire values. The accuracies 


of the experimental measurements of u/U and v,/U were estimated to be 


0 0 


about 0.5 percent and 1.0 percent, respectively. 

Figure 7 shows the good agreement between the theoretical and the 
"X¥" hot-film experimental axial velocity profiles, and even better agreement 
with the theory for the single element sensor. Wang and Boers” modified 
the eddy viscosity model, and this modeling plays an important role in 
computations of axial velocity profiles near the wall. The agreement be- 
tween theory and experiment is shown in Figure 7 to be very satisfactory 
near the wall. 

The agreement between theory and experiment for the radial velocity 
profiles is generally less satisfactory. This velocity component is small 
and more difficult to measure accurately. Nevertheless, the maximum 
difference between the computed and measured values of v_/U is less than 


0 
0.03. 


COMPARISON OF MEASURED AND COMPUTED INTEGRAL PARAMETERS 
x 
The distributions of (radial) displacement thickness ©. and (radial) 
boundary layer thickness 6. are given in Figure 8 along with the stern 


profile of the axisymmetric body. The theoretical boundary layer thickness 


30 


€L8°O “TE8°O “¥0L°O = 1/X ‘syq8ueT TeTXy [euUCTSUeUTpuON - e/ eAN8Ty 


ail 
An. X 
v0 70 00 
ce ae te ean 
4 — 10 — 
95460 = — (a) tn Xn 
v0 Z0 
co eu 
Oe on (V7) 
v0 
841M 10y 1UeWe;e ejBuls * 2 O 
wy oy xX = S 9'0 
£6 @” 
SLNIWIN3d x3 © 
ni 80 
V 


“Ore 
€28°0 = 5 (9) 


S.NIWI83d x3 n 


son 
241M 10y 1UeWe;/e B/buls xm 


© 
r—) 


Slelat sais! 


. wiyiy0Y |X, “== O 
OL Xn _ 


eal ‘0 
L€8'0 =~ (a) b 7 (oo 


SIN3IWId3d x3 


ane. 
vOL0 =v) 


¢ Apoqieijy jo tohey ATepunog urzeqs ssor0y 
uotjnqtizastqg AILOOTeA [eTpey pue [erxXy ueoW peanses, pue peqnduog - / osan3Ty 


0:0 


c0 


v'0 


90 


8'0 


OL 


Sal 


EQ630)  LG6s0r 60680) = 1) = 


ao 
euIM 104 1UeWe/e B/BuUIs a 
n 


wy oy XK on 
x 
n 


a 
by 


SLNIWIW3IdXF 


‘syqsuey] [TeTxy [TeUOTSUSUTpUON - q/ oan3Ty 


Vv 
Vv 
Vv 
WE 
» 


OV 
Ky, 
ay, 

V 


41m 104 
epi 
juawajea ajbuis To @SIAAT OU es 
n 
jueweye ejBuis * = 
" > 
wily OY XK. = OL 
n 


A 


‘ ‘n 
wy 2OY XK. 


n 
4n 


SLIN3IWIH3d x3 


CL SiNawivadxa 
ZLG6°0 = 4 (3) 


bane 
£606'0 = + (d) 


(penutquo)) ¢ s1n3Ty 


32 


THEORETICAL EXPERIMENTAL 


STERN 
PROFILE 


“0.70 0.75 0.80 0.85 0.90 0.95 1.00 1.05 1.10 1.15 1.20 


Figure 8 - Comparison of Measured and Computed Boundary Layer 
and Displacement Thickness Parameters 


is defined as the normal distance from the body where the nondimensionalized 
tangential velocity u,/U, has a value of 0.995. In this definition, uy is 
the dimensional tangential velocity and Wi is the dimensional local 
tangential velocity calculated by the potential flow program using the 
latest body-wake displacement model. 

The measured value of the boundary layer thickness 6 is a difficult 
value to obtain. This variable is defined as the radial position where 
Reynolds stress uv! becomes "small." The overall accuracy of the measured 
values of OF presented in the present investigation is estimated to be 
about 10 percent. 


The integral parameters of interest are: 


ae Csi Ua reee is 
aperin) Yad 2s r 
a i oleeny, ® 
max thle 6 max 
O° max 


Norn I ( =| ( r q ( r (9) 
Diy 7 o Ts r 
r 6 max max 
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(10) 


x 
— (11) 
Us 6 ee 


max 
if 


Te ifse 
O° max 


(rot) /T ax 


oe { (- 3 (=) EEN di 2 (12) 
r 4 7 as (a le) 


where Us is the value of uy Bye Te To 0 


value of boundary-layer thickness measured normal to the body surface. 


ate O. = fe, ap ns sec oO and ns is the 
The values of integral parameters derived from the measured axial velocity 
profiles, u,/Us, are given in Table 4. Figure 8 shows the theoretical and 
measured values of Cradiat) boundary-layer thickness 6. and the (radial) 
displacement thickness a The agreement between measured and computed 


x 
values of 6. and 6. is satisfactory. 


MEASURED TURBULENCE CHARACTERISTICS 
The cross-wire probe was used to measure the turbulence character- 
istics in the thick boundary layer. The measured Reynolds stresses and 
the measured mean velocity profiles were used to obtain eddy viscosity and 


mixing length. 


MEASURED REYNOLDS STRESSES 


The exchange of momentum in the turbulent mixing process is repre- 

. ° ° ' y 12 2 

sented by the distribution of the Reynolds stresses -u'v’, u' , v' , and 
We 92 WZ We ; ; ; 

W where u' , v' , and w are the turbulent fluctuations in the axial, 
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radial, and azimuthal directions, respectively. With the probe elements 
alined vertically, the Reynolds stress and the axial and radial fluc- 


tuations are computed by an on-line computer with a sampling rating of 


1024 data values in 8 seconds. To obtain the azimuthal fluctuation un 
the hot film sensor is rotated 90 degrees. Data points are accepted at the 
same rate and analyzed by the same computer. 

Figure 9 shows the nondimensionalized measured distributions of 


Eee We) F 
Reynolds stress -u'v'/U. and the three components of turbulent fluctuations 


0 


y u'?/U5, Vv" /Up» andy w'*/U, at various axial locations along the 


axisymmetric body. The axial fluctuation component is the largest of the 


three components for all axial locations. The radial velocity fluctuation 
component has the smallest value. The intensity of each of the three 
components is larger near the surface of the stern and reduces as the 
boundary-layer thickens for x/L < 0.95. A sharp gradient exists in the 
turbulence velocity components near the surface since velocities go to 
zero at the wall. 

The maximum value of the nondimensionalized Reynolds stress av" /U5 
occurs near the wall for x/L < 0.873 with a strong variation in the radial 
direction.* As the axial location increases and the boundary layer thick- 
ends, the radial location of maximum Reynolds stress moves away from the 
wall. The Reynolds stress reduces quickly from the maximum value to zero 
at the edge of the boundary layer. At axial locations of x/L = 0.951 and 
0.987, the measured values of -u'v' reduce from the maximum value to 
zero at the wall. These results are consistent with those of Huang et al. 
for the two afterbody shapes tested. 


il F A ; 
Huang et al. also have presented measured distributions of a turbu- 


lence structure parameter a, where a, = elec and at She 0o a GI a Ge 


1 Hf 


Thin boundary layer data show a, to have a constant value of approximately 


1 
0.15 between 0.056 and 0.86. The data of Huang et cil,” indicated that for 


a thick axisymmetric boundary layer, had a value of approximately 0.16 


oat 


*The spatial resolution of the "X' hot-film probe used may not be fine 
enough to measure the Reynolds stresses precisely near the wall of the thin 
boundary layer. 
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Figure 9 (Continued) 
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Figure 9b - Nondimensional Axial Lengths, x/L = 0.951, 0.987 


up to (r-r9)/5_. = 0.6, and that the value of ay decreased toward the edge 
of the boundary layer. The values of ay measured on the present axisym-— 
metric body are shown in Figure 10. The present results are similar to the 
earlier results of Huang et ails The values of ay decrease in the inner 
region at the location of x/L = 0.987. 

No attempt was made to remove the free-stream turbulence fluctuations 
from the measured components of turbulence velocities. The reduction in 
the value of a, near the edge of the boundary layer is in part caused by 
the larger contribution of the free-stream turbulence to ae than to -u'v" 


(Reference 1). 


Figure 10 - Measured Distributions of Turbulent Structure Parameter 
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EDDY VISCOSITY AND MIXING LENGTH 
The experimental values of eddy viscosity e€ and mixing length & 


across the thick stern boundary layer are derived from the measured values 


du 
of Reynolds stress -u'v' and the mean velocity gradient = using the 
Ei hy 5 5 163) 
definitions 
du, 
nq v = sae 
u'v € a (13) 
and 
9 du, du, 
— v v = aie 1 aes ote 
u'yv Q aE Tea (14) 


where the Reynolds stress -u'v' and the axial velocity u, are measured by 
the "X" hot-film probe. By fairing the measured velocity profile uy. with 
a 'spline' curve, the value of du, / 9. can be computed numerically. 4 

The experimentally-determined distributions of eddy viscosity e/U. ©. 


are given in Figure 11 for the axisymmetric body under consideration. The 
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Figure 11 - Measured Distributions of Eddy Viscosity for 
Afterbody 5 
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solid curve in this figure is the thin boundary layer eddy-viscosity 
model of Cebeci and ene The experimental values of eddy viscosity for 
axial locations, x/L < 0.873, where the boundary layer is thin, agree 
reasonably well with the Cebeci and Smith model. However, as the boundary 
layer thickens, the measured values of eE/U. become nearly one-fifth o 
the values for thin boundary layers given by the Cebeci and Smith model. 
The experimentally-determined distributions of mixing length 4/6 are 
shown in Figure 12. The solid curve in this figure gives the thin boundary 
layer results of Bradshaw, Ferriss, and ewes ue Again, agreement is 
reasonable for x/L < 0.873 and the experimental values in the thicker 
boundary layers are roughly one-third the values in the thin boundary 
layers. Similar reductions of eddy viscosity and mixing length also were 
measured earlier by Huang et ales and by Patel and Tyee)? The 
experimentally-derived eddy-viscosity and mixing-length data are also 


given in Table 4. 
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Figure 12 - Measured Distributions of Mixing Length for 
Afterbody 5 
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TURBULENCE MODEL 

Figure 12 shows that the distribution of measured mixing length in 
the thin boundary layer is roughly approximated by the thin boundary layer 
results of Bradshaw et Bae However, the thick boundary layer values of 
4/6 are considerably smaller than Bradshaw's thin boundary layer results. 
Figure 12 also shows that the maximum values of the mixing length occur at 
approximately the same radial location for theory and experiment. Huang 
et ale proposed that the mixing length of an axisymmetric turbulent 
boundary layer is proportional to the square root of the entire turbulence 


annulus between the body surface and the edge of the boundary layer 


y D2) 
Rox (r9t5_) “6 (15) 


The present measured values of 2/Y (rt6) =r 5 are given in Figure 13 
where the solid curve is derived by Huang et ails This figure shows that 


the values measured agree reasonably well with Huang's hypothesis. 


Figure 13 - Similarity Concept for Mixing Length of Thick 
Axisymmetric Stern Turbulent Boundary Layer 
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The new mixing length definition has been incorporated into the 
Douglas C-S boundary layer method of Wang and Fnemes° The original mixing 
length formulation is approximately valid and is used until the boundary 
layer thickness is greater than 20 percent of the body radius (670.2 ro): 
At this point, the boundary layer is considered thick and the mixing 


length 2 is defined by 


LS LS (16) 


Huang derived this equation and the limits on its use from the maximum 


values of mixing length given by Bacleberrs Los and the experimental data 


of Huang et aie 


TURBULENCE REYNOLDS NUMBER 


The axial turbulence velocity Reynolds number is defined as 


fur? 
aE @ Ce 


(17) 


where A is the microscale of the turbulence computed by the relationship 


2 = (18) 


Here, Taylor's hypothesis, the assumption of an isotropic field of turbu- 
lence (although not quite true in Figure 10) and space-time equivalence, 
is used to obtain the microscale X}. The time differentiation of ul was 
accomplished by using an analog operational amplifier (as used by Frenkiel 


et aia) in a differentiating mode. The departure from linearity as a 
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function of frequency was found to be 1 percent at 10 kHz and 5 percent at 
20 kHz. The axial velocity fluctuation ul was recorded on magnetic tape 
for 30 sec at a tape speed of 152.4 cm/s using a multichannel recorder. 
This tape recorder had a rapid fall-off in response above 20 kHz. The 


recorded data for ul were played back to determine the mean square values 


of 0 and (Surat) (with the signal passing through the differentiation 
x 
amplifier) by a true RMS voltmeter. 
The measured values of turbulence Reynolds number R) at various axial 


locations of Afterbody 5 are shown in Figure 14. The value of Ry decreases 


0 100 200 300 400 500 600 700 800 300 1000 


Figure 14 - Measured Turbulent Reynolds Numbers 


with increasing axial distance x, which is a distinguishing characteristic 
of a thick stern boundary layer. For a given x/L, the measured values of 
Ry are almost constant for distances less than 60 percent of the boundary 
layer thickness. The measured values of Ry in this thick San 
layer are higher than those generated by most turbulence grids. The 


turbulence Reynolds number Ry in the present thick stern boundary layer is 
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sufficiently high that the results obtained should be typical of high 


Reynolds number thick stern turbulent boundary layer flows. 


TURBULENT LENGTH SCALE IN THICK STERN BOUNDARY LAYER 
The experimentally derived mixing length data in the thick stern 
boundary layers of Afterbodies 1, 2, and 5 indicate that the mixing-length 
scale is proportional to the square root of the annular area between the 
body surface and the edge of the boundary layer. It is of fundamental 
interest to perform an independent and more direct measurement of the 
turbulence length scale to support this empirical result. The spatial 


correlation coefficient of two axial velocity fluctuations at points A and 


B, Pap (4) = aren) fat? Vay as a function of the radial separation 
distance between the points A and B, al Sea Laie is chosen for this 
purpose. The axial velocities at points A and B were measured by position- 
ing two single-element TSI (Model 1212) hot-wire probes with 90-deg bends 
at the two points. The linearized outputs of the velocities from the two- 
channel TSI (Model 1050-1) hot-wire anemometer were a.c. coupled by using 
an amplifier with low frequency rolloff set at 1 Hz. The fluctuating 
velocity correlation aon was then measured by a TSI Correlator (Model 
1015C). The spatial correlation coefficients were measured at various 
axial locations in the thick stern boundary layer of Afterbedy 5. At 


each axial location, point A was fixed at r, = r, + 0.167 6. and the radial 


A 0 
location (2,78 4) of point B was increased until the value of Pap became 
MySTiicelel ea 

The measured spatial correlation coefficients are plotted in Figure 15 


(data are tabulated in Table 5) against three nondimensional length 


parameters, namely, n = d_/¥ (rytd_)?-15 d/o» and a./8)- As shown in 
Figure 15, the measured correlation coefficients at various axial locations 
are not similar when n is taken to be d/o. or a,/833 the correlation 
length scale of the turbulence in the thick stern boundary layer is not 
proportional to the local radial boundary-layer thickness é. or the planar 


x 
displacement thickness O02 However, it is interesting to note that the 
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Figure 15 - Similarity Length Scale in Thick Axisymmetric Stern 
Boundary Layer 
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Figure 15a - Nondimensionalized by Turbulent Annular Area 


Figure 15b - Nondimensionalized by Planar Displacement Thickness 
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Figure 15 (Continued) 


Figure 15c - Nondimensionalized by Boundary Layer Thickness 


measured correlation coefficients at various axial locations in the thick 


stern boundary layer do possess a similarity form when the nondimensional 


length parameter n is taken to be al (eot6,) =r. This independent and 
more direct measurement of the turbulence length scale as determined from 
the spatial correlation coefficient measurements provides additional support 
that the appropriate length scale in the thick axisymmetric stern boundary 
layer is proportional to the square root of the turbulence annular area 
between the body surface and the edge of the boundary layer. Throughout 

the entire stern boundary layer, the measured Taylor microscale 4 at each 
axial location is found to be approximately equal to the measured corre- 
lation length scale ds for which the measured correlation coefficient 


Pap (d_) is equal to 0.4 (Figure 15). The measured values of 


A/ (ryt6,)°-r6 were found to vary from 0.024 to 0.036 for the five axial 


locations (Figure 15). This is additional evidence of the similarity 


46 


96T°0O 


TIT‘O 


09T‘O 909°O OvT"O Ss0°0 79T°O 
(A a0) 91S°0 €€T'O 8S0°0O 87T 0 
£€T "0 sos ‘0 99IT‘O 860°0 £€T°O 
SOT‘O 107 °0 9760°0 6TT°O cIT"O 
£960°0 99€°0 S780°0 9ST°O €S60°0 
0280°0 TEE*O 79L0°0 S6T*O €620°0 
7690°0 £92°0 4090°0 ¢lz°0 £990°0 
71S0°0 461T°O 6770°O €7£°0 89470°0 
6170°0 6ST°O 89£0°0 €17°0 0se0°O 
0920°0 £860°0O 8220°0 TSS °0 79720°0 
= 
d 05. (79+°3) A 
9 (G4 
¥ fv, 
as. Ey 


oss"0 
£6%7°0 
977 °0 
7LE°O 
0zE*0 
997°0 
t7z ‘0 
LST°O 
LTT"0 
7880 °0 


d 
g 


ples 
I 


Szt‘°O 
€TT ‘0 
STOT‘O 
1S80°0 
L£7L0°0 
s090°0 
90S0°0 
LS€0°0 
£920°0 
T0Z0°0O 


"610  6£6°0 0Sot‘0 GOO OES Nae SLET‘O 
o9T'0 = ELL" 7980°0 040°0 ZLz°0 ~—Yh6°T STIT‘0 
EYT'0 -€69°0 $£L0°0 460°0 O70 O049°T T%60°0 QOr-O} Slezo Teere 21°0 
ZIT‘0 ZS *0 9090°0 z9T"0 = ELT'0.-S «EZ 6020°0 LyT'0 Z0E°0 ~—-989°Z 00T'0 
$960°0  L9%°0 zzS0°0 THEO SO) POE 8190°0 O8T'O Shz"0 © 9LTZ 0180°0 
7080°0  68€°0 S€0°0 €0£°0  SOT°O 9h2°0 8740°0 z9z"0.-=-s STO ~—2T9°T 0090°0 
1”90°0 —ETEO 0S€0"0 0Sr°0 T090°0 ~—«6zv"0 9200 OSS Waa EOE 00%0°0 
40%0°0  L61°0 02z0°0 09S°0 40'0  4TE"0 O8TC *0 oSr'0 ~=— #L0"0-~Ss«S 910 620°0 
5 Co 92-9492) A = d 92-7 (943) gS 9a-e+a) A 
a) ee Sane eee ee || EUW at ee a, pee Sees 
1 a ef i 1a 1a a 1 
P P P P P P P P 
-_—— 
eur aa, aaa 23, qe ‘ Zen xeu, xeu 
26050) = —7— S970) = = WO OS OS 0°0 = ‘Onz"0 = 
9 9 9 9 as eee = i 
* * x? 9 
eae Nie are eer Airs [E7g@ oH ae 1 
Bes"0 = > 606°0 == 86L°0 = 80-4 916°0 = — ‘tes0 = 2 


SNOTLONN4 NOILVIHYYOO IVILVdS CHYNSVAW - 


G ATEVL 


47 


property of the turbulence in the thick axisymmetric stern boundary layer. 
The present set of experimental results is consistent with the conjecture 
of Townsend that the flow patterns in a wake and in a boundary layer 
possess an overall structural similarity. It may not be necessary to use 
"local" equations for turbulent Reynolds stresses with a diffusion term 


aa 19. 
and with a correction to the extra rate of strain. 


CONCLUSIONS 

The present investigation of the thick axisymmetric turbulent 
boundary layer on an inflected stern without shoulder separation is an 
extension of the earlier research reported by Huang et ail," An 
additional comprehensive set of boundary layer measurements, including 
mean and turbulence velocity profiles, turbulence Reynolds stresses, static 
pressure distribution, and correlation length scales of the turbulence, are 
presented. The following major conclusions can be drawn. 

The Beeston! and Tile enna displacement body concept has further been 
proven experimentally to be an efficient and accurate approach for computing 
the viscid and inviscid stern flow interaction on axisymmetric bodies. 

The measured static pressure distributions on the body and across the 
entire thick boundary layer were predicted by the displacement body method 
to an accuracy within one percent of dynamic pressure. 

Neither the measured values of eddy viscosity nor the mixing length 
were found to be proportional to the local displacement thickness or the 
local boundary-layer thickness of the thick axisymmetric boundary layer. 
The measured mixing length of the thick axisymmetric stern boundary layer 
was again found to be proportional to the square root of the area of the 
turbulent annulus between the body surface and the edge of the boundary 
layer. This simple similarity hypothesis for the mixing length and the 
displacement body concept have been incorporated into the Douglas C-S 
differential boundary-layer method? by Wang and Euareme The improved 
method predicts satisfactorily the measured mean velocity distributions in 
the entire stern boundary layer. 

The measured correlation length scales of the turbulence and the 


measured Taylor microscales in the entire stern boundary layer are found to 
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be proportional to the square root of the annular area between the body 
surface and the edge of the turbulent boundary layer, which indeed is the 


appropriate overall length scale of the thick stern boundary layer. 
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